Inversion of Optical Reflectance in the Fullerenes 



F. Marsiglio 

Department of Physics, University of Alberta, Edmonton, AB T6G 2J1 

Since the discovery of superconductivity in the alkali-doped fullerenes JlJ, the electron-phonon 
interaction has been the primary suspect for superconductivity in this class of compounds. In this 
paper we first provide a pedagogical review of how the question of mechanism has traditionally been 
settled, and then some well-known properties of the optical reflectance (and the derived conductance) 
are summarized. Finally we demonstrate how a recently derived inversion procedure can use the 
optical properties of a metal to infer the magnitude of the electron-phonon interaction. We conclude 
that in the case of the alkali-doped fullerenes, the electron-phonon interaction is sufficiently strong 
to explain the transition temperatues found in these the materials. 



Introduction 

The fee solid composed of Cgo molecules at the lat- 
tice sites is an insulator, and almost certainly a band 
insulator ||. As is the case with graphite, this mate- 
rial can be doped with alkali metal atoms. One can try 
for various 'stoichiometric' results, such as A„C6o, where 
A is an alkali metal atom, and n is an integer ranging 
from 1 through 6. Various interesting phases result, as a 
function of n (and to a lesser extent A). In some cases, 
the compound apparently doesn't form, with indications 
of phase separation at low temperature (see Ref. || , for 
example). The point of interest for us is n = 3, for it 
is at this doping level (and apparently only this doping 
level) that superconductivity occurs. Unlike the focus of 
much of the work in high temperature superconducting 
cuprate materials (which also have an insulating parent 
state), it appears to be a much more difficult problem 
here to relate the various phases that occur as a func- 
tion of n. Here the situation is intrinsically different, in 
that only stoichiometric compounds form (although sug- 
gestions that this may also be the case in the cuprates 
have appeared from time to time). Thus, for this work 
we relinquish the 'big picture', and focus only on A3C60, 
and in particular, A = K, which is a superconductor with 
a transition temperature, T c = 20 K. 

There are many indications |3j^|] that the K3C60 com- 
pounds are (a) well-behaved Fermi Liquids, and (b) con- 
ventional electron-phonon superconductors. Property (a) 
does not imply property (b), although property (b) pre- 
supposes that (a) is true. The main indication that 
Fermi Liquid theory (FLT) might be the appropriate 
starting point is the temperature dependence of the re- 
sistivity, which exhibits the standard FLT behaviour, 
p{T) — a + bT 2 M. Various measurements below the crit- 
ical temperature indicate that the superconducting state 
is a conventional s-wave one. Gap measurements, using 
both tunneling and infrared spectroscopy ^-^| indicate 
an s-wave gap, though there is some controversy over the 
size of the gap. The existence of an s-wave gap is also 
supported by penetration depth measurements ]To| . Fur- 
thermore, a coherence peak in the spin relaxation rate, 



1 /T\ , has been measured using muons Q , a result which 
lends strong support to an s-wave BCS picture. Finally, 
the isotope effect has been determined to be significant 
(see Ref. for a summary), although again the size of 
the isotope exponent is controversial. 

The Question of Mechanism 

All of these measurements support the existence of an 
s-wave gap, but, with the exception of the isotope ef- 
fect, have little to say about mechanism. Partly this is 
the "curse" of universality; BCS theory was so successful 
partly because it provided universal results (for various 
properties and temperature dependences). However, one 
consequence of this is that for a wide variety of mecha- 
nisms, if the coupling is weak, the prediction for many su- 
perconducting properties is the same. Hence, the theory 
is not very useful for discriminating mechanism. One can 
claim, for example, that in the case of Al (which is a very 
weak coupling superconductor) , there is no definite proof 
for electron-phonon superconductivity E3j. For certain 
superconducting materials (known in the older literature 
as "bad actors") like Pb and Hg, deviations from uni- 
versal behaviour were measured, so there was hope that 
these deviations could point towards a mechanism. This 
is exactly what happened in the case of Pb, particularly 
through tunneling experiments, as we will now describe. 

In single particle tunneling experiments an (insulating) 
oxide layer is sandwiched between a superconducting ma- 
terial and a normal metal (other variants are possible, 
but we focus on this, the simplest). An external circuit 
allows electrons to tunnel across this barrier. If a single 
particle energy gap exists in the superconducting side, 
then no current will flow until the bias voltage exceeds 
this energy gap. In this way one can measure the en- 
ergy gap in a superconductor (provided one exists). This 
is exactly what the original experiments by Giaever |l5| ] 
measured (in an Al/Pb sandwich). 

Actually, even this measurement indicated that some- 
thing was wrong. The gap ratio, 2A/kBT c , was measured 
to be well over 4 (later more accurate measurements indi- 
cated a value closer to 4.5), whereas BCS theory predicts 
a universal result, 3.5. Strong coupling was suspected, 
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although strong coupling extensions of BCS theory only 
managed to raise this universal value to yet another uni- 
versal maximum value, 4.0 |H. By this time Eliashberg 
had extended the theory to include dynamical interac- 
tions with phonons ]l5| |. As further refinements in both 
theory and experiment developed, it became clear that 
certain structure in the density of states (as measured 
by electron tunneling in Pb, one of the "bad actors") 
was due to the electron- phonon interaction JL6 17 1. Fi- 
nally, McMillan and Rowell succeeded in using tunneling 
measurements and the Eliashberg equations to extract 
the underlying electron-phonon interaction in the form 
of the electron-phonon spectral function, a 2 F(Q) |l8| , |l9| . 
We now describe qualitatively the essence of Eliashberg 
theory (for various reviews see Ref. [^p|-^3|), and of the 
McMillan-Rowell inversion procedure. 

Within BCS theory, the gap parameter AgQg serves 
both as an order parameter, and as the energy below 
which the density of states is zero (hence the name "gap" 
parameter). The only input required in the theory is 
a model attractive potential, multiplied by the electron 
density of states at the Fermi level: N(cf)V. As al- 
ready mentioned above, Eliashberg theory goes one step 
further; it recognizes that a source for the attraction is 
the electron-phonon interaction, and that retardation ef- 
fects can be important since the mean ion velocity is 
significantly less than the Fermi velocity of the electrons. 
The result is that the source of attraction is encapsulated 
by an electron-phonon spectral function, a 2 F(£l). More- 
over, the order parameter in Eliashberg theory is fre- 
quency dependent. Finally, the direct electron-electron 
Coulomb repulsion is recognized as a source of depairing, 
and is generally included at a minimal level through a 
single dimensionless parameter, fi* . The '*' is present 
to remind us that this is a pseudopotential, and ought to 
have a greatly reduced value if one uses cutoffs in the the- 
ory which apply to the phonon energy scale. To summa- 
rize, Eliashberg theory produces an equation analagous 
to the BCS one, except that the order parameter is a 
function of frequency u>, a functional of a 2 F(n), and a 
function of /i*. It is also complex (not in the usual sense: 
this applies even if the phase has been set to zero) . Once 
this gap function, A(w), is obtained, one can calculate 
a variety of properties, one of which is the single par- 
ticle density of states, N(uj) = ^V(ef)Rc 



This expression bears a remarkable resemblance to the 
weak coupling result, for which the gap function would 
simply be replaced by the gap parameter, Ag^g. Note 
that the I-V characteristic is proportional to the single 
particle density of states: dI(V)/dV oc N(V), so that an 
I-V measurement immediately yields ths single particle 
density of states. 

The procedure formulated by McMillan and Rowell is 
the following. One first measures the structure above 
the gap edge in dl/dV accurately. Next, one "guesses" 



an a 2 F(fl), and using this guess computes A(w), and 
hence dl{w)/ dV, using Eliashberg theory. The calculated 
function is compared to the experimental result, and cor- 
rections to the initial guess are computed. Corrections 
can be obtained through functional derivatives; the pro- 
cedure amounts to a Newton-Raphson method for func- 
tions (rather than for single variables). The new a 2 F(fl) 
is used through the procedure all over again, until the 
result converges. One then claims to have "measured" 
the underlying a 2 F{Q). That this procedure works is a 
priori not obvious. There seems to be no reason why the 
solution should be unique, but apparently most of the 
time it is. That the result is physically meaningful is ap- 
parent because of several tests. First, that the function 
comes out positive is already a triumph. Negative spec- 
tral functions are not ruled out in the iterative process, 
so the positivity indicates a physical result (the spec- 
tral function must be positive (for positive frequency)). 
Although not apparent from the above description, the 
inversion procedure requires knowledge of the density of 
states in the phonon region only. Nonetheless, structure 
occurs well beyond the phonon region. One can then use 
the spectral function with Eliashberg theory to predict 
the density of states at higher energies. The agreement 
with the experimental result is spectacular (see Fig. 32 
in Ref. [l8[]). Finally one can compare to neutron scat- 
tering experiments, which measure the phonon density 
of states, F(Q), directly; these experiments tell us at 
what energies the phonons exist, but not how strongly 
coupled they are to the electrons. Hence, the two func- 
tions, F(Cl), and a 2 F(£l), don't have to be similar in 
many respects, though the frequency range of the latter 
clearly cannot exceed that of the former. As it turns out, 
these functions are rather similar (and, of course, fulfill 
the required condition); this indicates that the coupling 
function a 2 (fl) = a 2 F '(fi) / F '(f2) , is not very frequency 
dependent. A detailed comparison of these two probes is 
provided in Ref. Q . 

Finally we should mention that other measurements of- 
ten corroborate this story. The gap ratio often comes out 
in better agreement with experiment, and various ther- 
modynamic quantities show better agreement too |p3| . 
One might well wonder at how successful this approach 
has been historically, though, since important effects 
apparently have been omitted. For example, electron- 
electron Coulomb interactions have only been included 
insofar as they alter the underlying band structure (this 
is implicit), and through the single parameter /i* (this is 
explicit in the theory). 

Perhaps partially for the reasons mentioned in the pre- 
vious paragraph, along with other aspects, such tunnel- 
ing measurements have so far failed to provide conclusive 
results for the cuprate superconductors. In the case of 
the alkali-doped fullerenes, tunneling has clearly demon- 
strated a gap, but so far tunneling measurements have 
not been used to invert for a 2 F(fl) p5j. Part of the prob- 
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lem is that the phonons are known from neutron scatter- 
ing |^6| to extend out to about 250 meV, much further 
than in conventional superconductors. For these reasons 
we looked for an alternate method to deduce the under- 
lying interactions, one which uses the measured optical 
conductivity |^7j . In the next section we briefly describe 
optical measurements and then describe our procedure 
and its application to K^Ceo- 

Optical Properties in Metals 

A common method to determine the optical conduc- 
tivity is to measure the reflectance |Q as a function of 
frequency, usually at normal incidence. The reflectance, 
R{v), is defined as the absolute ratio squared of reflected 
over incident electromagnetic wave amplitude. The com- 
plex reflectivity is then defined by 



r{v) = R(v)exp (i6(v)), 



(1) 



where 9{v) is the phase, and is obtained through a 
Kramers-Kronig relation from the reflectance p3] 
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The complex reflectivity can be related to the complex 
index of refraction, n(u), 
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(3) 



which, finally, is related to the complex conductivity, 
<j(v) (using the dielectric function, e{v): 
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FIG. 1. The reflectance R(y) vs. frequency for the Drude 
model. Three curves are shown, corresponding to the three 
elastic scattering rates shown. 



where eoo is the dielectric function at high frequency (in 
principle, for infinite frequency this would be unity). It 
is through such transformations that the 'data' is often 
presented in 'raw' form. Nonetheless, assumptions are re- 
quired to proceed through these steps; for example, Eq. 
(||) indicates quite clearly that the reflectance is required 
over all positive frequencies. Thus extrapolation proce- 
dures are required at low and high frequencies; a more 
thorough discussion can be found in |29) ; see also |j(J . 

In Figs. 1 and 2 we illustrate how elastic and inelastic 
scattering processes affect the reflectance. In Fig. I we 
use a simple Drude form for the complex conductivity, 
o(y) = oo/(l — vt) and plot the reflectance which results 
from the application of the above formulae. Here the 
coefficient o~q is the dc conductivity, and the parameter 
1/r is the electron elastic scattering rate. The source of 
clastic scattering is unspecified, but usually results from 
electron-impurity scattering. In Fig. 1 the plasma edge 
is evident, and occurs, of course, near the plasma fre- 
quency, u>p. The "sharpness" of the dropoff decreases 
with 
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FIG. 2. The reflectance R(y) vs. frequency for the case 
shown in Fig. 1, but with additional inelastic scattering 
(modelled through the electron- phonon interaction in Pb). 
Note the general rounding of the curves and the decreased 
reflectance (as compared with those in Fig. 1). 
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increasing scattering rate. The plateau region below the 
plasma frequency is referred to as the relaxation region; 
the level decreases with increasing scattering rate, as 
expected, since stronger scattering should cause the re- 
flectance to decrease. Finally, at very low frequency, the 

Hagen- Rubens form, R(v) w 1 — ,/ — — — takes the re- 
flectance smoothly to unity at zero frequency. 

In Fig. 2 we use the same elastic scattering rates for 
the three curves shown in Fig. 1. In addition we include 
inelastic scattering through the electron-phonon interac- 
tion (in this case with Pb). Note that the main effect 
is a rounding of the edge, and a general decrease of the 
reflectance everywhere. The effect is very nonlocal in fre- 
quency — note that the Pb phonon spectrum extends to 
about 10 meV, which is tiny on the scale of this figure, 
and particularly on the frequency scale of the changes 
with respect to Fig. 1 . The effect of the inelastic scatter- 
ing is also cumulative, a result which follows from the fact 
that the effective scattering rate increases as a function 
of frequency J3l],|3(| . 

The corresponding conductivity curves are shown in 
Figs 3-5. In Fig. 3 we show the Drude form (Lorentzian) 
for the three scattering rates shown previously. The 
results shown in Fig. 3 are independent of tempera- 
ture. In Fig. 4 we suppress the elastic scattering rate 
(1/r = 0), to see more clearly the effect of inelastic scat- 
tering (here included through the electron-phonon inter- 
action in Bai-zK^BiOs, with T c — 29 K). The results 
shown are in the normal state only; we have referenced 
the temperatures to the critical temperature as a conve- 
nience only. What is clear is that at high temperatures, 
enough phonons are present to mimic impurity scatter- 
ing, so that a "Drude-like" peak is present. At low tem- 
peratures this peak vanishes, and structure becomes more 
visible (solid curve in the figure). This structure provides 
a signature of the electron-phonon interaction, analogous 
to the information inherent in the structure in the tun- 
neling I-V characteristic [ ^2|j3^ ]. Finally, in Fig. 5, we 
show similar results, but with an elastic scattering rate 
1/t = 10 meV, so that a Drude peak remains even at 
zero temperature. Note that the structure remains, but 
is not nearly as visible as in the clean limit, and neither 
seems to display the degree of structure seen in tunnel- 
ing. We also remark that similar statements apply to the 
imaginary part of the conductivity (see Ref. Q for a full 
discussion) . 
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FIG. 3. The real part of the optical conductivity in the 
Drude model (elastic scattering only) for the same three scat- 
tering rates shown in Fig. 1. These results are independent 
of temperature. 
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The Inversion of the Optical Conductivity 

The inference of the electron-phonon interaction 
through the optical conductivity has a long history. No 
doubt it was apparent that the information was avail- 
able right from the initial studies of Holstcin |35| , where 
he made clear that inclusion of the electron-phonon in- 
teraction would have a significant impact on the optical 
conductivity. The first real attempt, however, to infer 



FIG. 4. The real part of the optical conductivity in the 
clean limit (no elastic scattering). The only source of scat- 
tering is inelastic, as modelled by the electron-phonon inter- 
action in Bai-zKzBiOs (T c = 29 K) The curves shown 
are all for the normal state. Note the "Drude-like" peak at 
the origin at high temperatures, and its absence as the tem- 
perature approaches zero. 
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FIG. 5. The real part of the optical conductivity as in Fig. 
5 but now with additional elastic scattering (1/t = 10 meV). 
A Drude peak is present at the origin for all temperatures. In 
both this figure and in Fig. 4 structure is only discernable at 
the lowest temperatures. 

electron-phonon effects from the optical conductivity was 
made by Allen J3(|. He was motivated by infrared mea- 
surements on superconducting Pb [j37| , and the possibil- 
ity of inversion in that case. The result was a qualitative 
success, and even yielded quantitatively good results |38| . 
Our recent work, which we now describe, builds on the 
work of Ref. §(J. 

The remainder of this section concerns only the normal 
state, in contrast to tunneling, where the superconduct- 
ing state was a requirement. This is both an advantage 
and disadvantage, as the theory we present is at zero tem- 
perature, which may be difficult to achieve in the normal 
state for some materials. The required expressions are 
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is the effective electron self-energy due to the electron- 
phonon interaction. The spectral function that appears 
in Eq. (^) is really a closely related function, as has been 
discussed by Allen j3(| and Scher |39). For our purposes 
we will treat them identically. Eqs. (||) and eluci- 
date why the electron-phonon signature is rather modest 



in the optical conductivity (it is even better hidden in 
the reflectance); two integrations effectively average over 
the electron-phonon spectral function. One would like 
to "unravel" this information as much as possible before 
attempting an inversion, so that, in effect, the signal is 
"enhanced" . To this end one can attempt various manip- 
ulations. As a first step one can make a weak coupling 
type of approximation to obtain |27fl the explicit result: 
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Insofar as Eq. (0) works extremely well (as we shall 
see in a moment), it is a remarkable result. It tells us 
that, with a judicious manipulation of the conductivity 
data, the underlying electron-phonon spectral function 
emerges in closed form. This is a much better result than 
we were initially hoping for. Not surprisingly it requires 
two derivatives of the data (recall the two integrations), 
so in practice accurate measurements are required (ac- 
tually smoothing the data works quite well — and too 
much smoothing cannot really give you a spurious result; 
at worst it will obscure an otherwise noisy result — see 
Ref. for an early application of Eq. (M)). Perhaps the 
most inaccurate part is the determination of the plasma 
frequency, which can be obtained through a sum rule or 
from other measurements, procedures which are both of- 
ten fraught with errors. 

Before we evaluate Eq. (Q) with a known example, 
we note that a full numerical inversion of Eqs. (j^) and 
(^) is possible. However, a straightforward application 
of a Newton-Raphson iteration technique turns out to be 
very unstable. We derived instead the following expres- 
sion Q: 
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so that the right hand side depends (in many places 
through both \(u>) and E(w)) on a 2 F(il). We have found 
that only a few iterations of these equations are needed, 
even with an initial guess which is blank. The result from 
this procedure in the case of Pb is plotted in Fig. 6 
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FIG. 6. a 2 F{v) for Pb (solid curve) 



along with the 



estimates obtained from Eq. (|7j) with an impurity scattering 
rate, 1/r = 1 meV (dotted) and 10 meV (dot-dashed). These 
are both qualitatively quite accurate, before they become neg- 
ative at higher frequencies. Also plotted is the result (dashed 
curve, indiscernible from the solid curve) obtained from a full 
numerical inversion, as described in the text. 



Q (dotted curve), where we have taken the coupling 
strengths and frequencies of particular phonon branches 
analyzed there, and arbitrarily broadened them. 

Our results are in very good qualitative agreement with 
the neutron scattering data |2(| . The energy scale is cer- 
tainly correct, and moreover the peaks line up fairly well, 
suggesting that the coupling is not dependent on energy. 
Note that we have omitted some negative pieces, as these 
are expected to arise from the use of Eq. (^) (not the 
low frequency negative component, however, for which 
we have no satisfactory explanation at present). We also 
note that the agreement with the photoemission-derived 
results is not particularly good. 

We can easily 'check' this result by computing the re- 
flectance expected from this spectral function and com- 
paring to the measured reflectance. The agreement is 
within experimental error p7[ . Thus, we have demon- 
strated a consistent explanation of the data (but not 
necessarily a unique one). The fitted parameters are 
u>p = 6000 cm -1 and 1/r = 95 meV for the plasma 
frequency and impurity scattering rate, respectively. 



(dashed curve), though the curve is essentially indis- 
tinguishable from the tunneling density of states (solid 
curve) which was used to generate the optical conductiv- 
ity data in the first place. This illustrates the success of 
the method. However, the more remarkable result is that 
of Eq. (Q). We have also plotted the result of this for- 
mula applied to the conductivity of Pb (again generated 
numerically) for two values of the elastic scattering rate 
(after all in a real experiment we don't know what this 
value is). The resulting curves (dotted and dot-dashed) 
give very good qualitative results (and don't depend very 
much on the value of 1 /r used) . They both exhibit a spu- 
rious negative part plus some spurious "recovery" noise, 
so these effects teach us about the limitations of Eq. (0) . 
If we ignore such negative pieces, we get a very good 
representation of the spectral function. 

The Case of Ks C§q 

For the case of K3C60 we simply use the approximate 
Eq. (Q) , bearing in mind some of the limitations already 
discussed. We have used the reflectance data at 25 K for 
K3C60 taken by Degiorgi et al. ||. Details of our analysis 
of this data are given in |^7j . The result of using Eq. (Q) 
is plotted in Fig. 7 (solid curve), along with the neutron 
scattering data from |^6| (dashed curve) . We also include 
a result obtained through analysis of photoemission data 
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FIG. 7. The a 2 F{v) for K 3 C 6 o (solid curve) extracted from 
the reflectance data of Degiorgi et al. js), using Eq. (Q). For 
purposes of analysis we have omitted the negative parts. The 
neutron scattering results from Ref. ]2^] (dashed curve) are 
also shown. Clearly the energy scale in o?F(v) matches that 
of the phonons, and some of the peaks even line up correctly. 
Finally, the dotted curve comes from an analysis of photoe- 
mission data p^ ] ; where we have arbitrarily broadened the 
phonon spectrum with Lorentzian lineshapes. 

We have used the extracted a 2 F(v) in Fig. 2 to deter- 
mine whether or not such a spectral function can account 
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for the superconducting properties in K 3 C6o- We find, 
given T c = 19 K, the Coulomb repulsion pseudopotential 
is fi* = 0.4 (using a cut-off, w c = 1 eV). This is in effect 
another check, since a negative /x* would indicate a break- 
down. The spectrum is characterized by A = 1.2 and 
u)i n = 40 meV, parameters which give BCS-like results 
for superconducting properties (for example, the gap ra- 
tio comes out to be 2Ao/1cbT c ss 4). If the low frequency 
peak in a 2 F(v) at 5 meV is excluded, the results are 
then very BCS-like: A = 0.8, [i* = 0.34, T c /uj ln » 0.016, 
and the gap ratio is 2Ao/ksT c « 3.6. These results are 
consistent with some tunneling and infrared measure- 
ments but inconsistent with others 0. The degree 
of coupling roughly agrees with that extracted through 
an analysis of photoemission data 0]. Microwave mea- 
surements may provide a more discriminating probe for 
the low frequency part of the spectrum, as discussed in 

In summary, we have described an inversion scheme to 
extract from the optical conductivity (or reflectance) a 
boson spectral function responsible for inelastic scatter- 
ing in a metal. The full inversion scheme is numerical, 
but we have also derived an explicit expression which pro- 
vides a remarkably good semi-quantitative result. We 
have applied this technique to K3C60, and have found 
the qualitatively correct a 2 F(v) for this material. Its 
frequency range lies in the phonon region and peaks in 
a 2 F(v) line up with peaks in the phonon distribution 
function as determined by neutron scattering. The cou- 
pling is sufficiently strong to explain the superconductiv- 
ity in this material. Thus, the weak coupling approach 
presented here appears to be able to account for the main 
features observed in the far-infrared. K3C60 appears to 
be a conventional weak coupling electron-phonon super- 
conductor. 
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